Dissociation of vertical semiconductor diatomic artificial molecules 
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We investigate the dissociation of few-electron circular ver- 
tical semiconductor double quantum dot artificial molecules 
at T as a function of interdot distance. Slight mismatch 
introduced in the fabrication of the artificial molecules from 
nominally identical constituent quantum wells induces local- 
ization by offsetting the energy levels in the quantum dots by 
up to 2 meV, and this plays a crucial role in the appearance of 
the addition energy spectra as a function of coupling strength 
particularly in the weak coupling limit. 
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Semiconductor quantum dots (QD's) are widely con- 
sidered as artificial atoms, and are uniquely suited 
to study fundamental electron-electron interactions and 
quantum effects There are many analogies with 'nat- 
ural' atoms. One of the most appealing is the capabil- 
ity of forming molecules. Indeed, systems composed of 
two QD's, 'artificial' quantum molecules (QM's), coupled 
either laterally or vertically, have recently been investi- 
gated experimentally ||,^ and theoretically Nev- 
ertheless, the direct observation of a systematic change 
in the addition energy spectra for few-electron (number 
of electrons, N < 13) QM's as a function of interdot 
coupling, has not been reported, and calculations of QM 
properties widely assume a priori that the constituent 
QD's are identical Special transistors incorporat- 

ing QM's 1^ made by vertically coupling two well defined 
and highly symmetry QD's Q are ideally suited to ob- 
serve the former and test the latter. 

In this work we present experimental and theoretical 
addition energy spectra characterizing the dissociation 
of slightly asymmetric vertical diatomic QM's on going 
from the strong to the weak coupling limits that corre- 
spond to small and large interdot distances, 6, respec- 
tively. We also show that spectra calculated for sym- 
metric diatomic QM's only resemble those actually ob- 
served when the coupling is strong. The interpretation 
of our experimental results is based on the application 
of local-spin density-functional theory (LSDFT) 0-^2). 
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It follows the development of the method thoroughly de- 
scribed in Ref. |12|, which includes finite thickness effects 
of the dots, and uses a relaxation method to solve the 
partial differential equations arising from a high order 
discretization of the Kohn-Sham equations on a spatial 
mesh in cylindrical coordinates |p^ . Axial symmetry is 
imposed, and the exchange-correlation energy has been 
taken from Perdew and Zunger |ll0| . 

The molecules we study are formed by coupling, quan- 
tum mechanically and electrostatically, two QD's which 
individually can show clear atomic- like features |^,^ . For 
the materials we typically use, the energy splitting be- 
tween the bonding and anti-bonding sets of single parti- 
cle (sp) molecular states, AsaSj can be varied from about 
3.5 meV for b = 2.5 nm (strong coupling) to about 0.1 
meV for b — 7.5 nm (weak coupling) [||. This is expected 
to have a dramatic effect on the electronic properties of 
QM's Figure 1 shows (a) a schematic diagram of 

a sub- micron circular mesa, diameter D, containing two 
vertically coupled QD's, and (b) a scanning electron mi- 
crograph of a typical mesa after gate metal deposition. 
The starting material, a special triple barrier resonant 
tunneling structure, and the processing recipe are de- 
scribed elsewhere HJlJl- Current /d flows through the 
two QD's, separated by the central barrier of thickness b, 
between the substrate contact and grounded top contact 
in response to voltage Vd apphed to the substrate, and 
gate voltage Vg. The structures are cooled to about 100 
mK and no magnetic field is applied. 

To analyze the experiments we have modeled the QM 
by two axially symmetric QD's. The QM is confined 
in the radial direction by a harmonic oscillator poten- 
tial muP'r^ jl of strength hu! — 5 meV (a realistic lateral 
confinement energy for a single QD in the few electron 
limit ||,|l^), and in the axial (z) direction by a double 
quantum well structure whose wells are of same width w, 
and have depths Vq ± S, with (5 < Vo Figure |l| (c) 

schematically shows the double quantum well structure 
and its unperturbed bonding and anti-bonding sp wave- 
fimctions. We have taken Vq = 225 meV and w = 12 
nm, which are appropriate for the actual experimental 
devices. If S is set to zero, the artificial molecule is 
symmetric ('homonuclear' diatomic QM); otherwise, it 
is asymmetric ('heteronuclear' diatomic QM). In the cal- 
culations here, (5 is 0, or is set to a realistic value of 0.5 
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or 1 meV In the homonuclear case AgAs is well re- 
produced by the law Asas(&) = Aq e~^/^° with 60 = 1.68 
nm, and Aq — 19.1 meV. It is easy to check that in 
the weak coupling limit 25 is approximately the energy 
splitting between the bonding and anti-bonding sp states 
which would be almost degenerate if 6 is 0. For this rea- 
son we call the mismatch (offset) the quantity 25. 

Figure ^(a) shows calculated addition energy spectra, 
A2{N) = U{N+l)-2U{N) + U{N-l), for homonuclear 
QM's with realistic values of h conveniently normalized 
as A2(iV)/A2(2). U{N) is the total energy of the A^- 
electron system. A2(A'^) can reveal a wealth of informa- 
tion about the energy required to place an extra electron 
into a QD or QM system |,|l|l. For small b (Asas ^ Tt^) 
the spectrum of a few electron QM is rather similar to a 
single QD, at least for TV < 7 At intermediate dot 
separation, the spectral pattern becomes more complex. 
However, a simple picture emerges at larger interdot dis- 
tances when the molecule is about to dissociate. For 
example, at 6 = 7.2 nm strong peaks at = 2, 4, 12, 
and a weaker peak at = 8 appear that can be easily 
interpreted from the peaks appearing in the single QD 
spectrum. The peaks at iV = 4 and 12 in the QM are 
a consequence of symmetric dissociation into two closed 
shell (magic) N = 2 and 6 QD's respectively, whereas 
the peak at iV = 8 corresponds to the dissociation of the 
QM into two identical stable QD's holding four electrons 
each filled according to Hund's first rule to give maxi- 
mal spin [^H]. The QM peak at iV = 2 is related to the 
localization of one electron on each constituent dot, the 
two-electron state being a spin-singlet QM configuration. 

Since the modeled QM is homonuclear, each sp wave- 
function is shared 50% — 50% between the two con- 
stituent QD's. Electrons are completely delocalized in 
the strong coupling limit. As b increases, Asas decreases 
and eventually bonding, |S), and anti-bonding, |AS), sp 
molecular states become quasi-degenerate. Electron lo- 
calization can thus be achieved combining these states as 
(|S)±|AS))/V2. 

We conclude from Fig. g(a) that the fingerprint of 
a dissociating few-electron homonuclear diatomic QM is 
the appearance of peaks in A2(A^) at TV = 2, 4, 8 and 12 
. This is a robust statement, as it stems from the well 
understood shell structure of a single QD. If we now com- 
pare this picture with the experimental spectra shown in 
Fig. ||(b), we are led to conclude that the experimental 
devices are not homonuclear, but heteronuclear QM's. 

The origin of the mismatch is the difficulty in fabri- 
cating two perfectly identical constituent QD's in the 
QM's discussed here, even though all the starting ma- 
terials incorporate two nominally identical quantum 
wells. This mismatch can clearly influence the degree 
of delocalization-localization, and the consequences will 
depend on how big 25 is in relation to Asas ■ Else- 
where we will discuss how the effective value of Asas is 
measured, and the mismatch is determined for all values 



of b |T^, so we merely note here that 25 is typically 0.5 
to 2 meV and nearly always with the upper QD (near- 
est top contact of mesa) states at higher energy than the 
corresponding lower QD states [see Fig. [^(c)]. 

Figure ^(b) shows experimental spectra, also normal- 
ized as A2(A^)/A2(2), for QM's with b between 2.5 and 
7.5 nm, deduced accurately from peak spacings between 
Coulomb oscillations (Id — Vg) measured by applying 
an arbitrarily small bias (Vd < 100 /iV). Likewise, also 
shown is a reference spectrum for a single QD [Q. The 
diameters of the mesas lie in the range of 0.5 to 0.6 /Ltm, 
and while all mesas are circular, we can not exclude the 
possibility that the QM's and QD's inside the mesas may 
actually be slightly non-circular, and that the confining 
potential is not perfectly parabolic as N increases [|l8|,|o| . 
We emphasize the following: i) The spectrum for the 
most strongly coupled QM {b ^ 2.5 nm) resembles that 
of the QD up to the third shell {N = 12). ii) For interme- 
diate coupling (b = 3.2 to 4.7 nm), the QM spectra are 
quite different from the QD spectrum, and a fairly no- 
ticeable peak appears at = 8. iii) For weaker coupling 
(b = 6.0 and 7.5 nm) the spectra are different again, with 
prominent peaks at N = 1 and 3. 

We confirmed the heteronuclear character of the QM's 
by performing LSDFT calculations with a 2 meV mis- 
match. The results are displayed in Fig. ||(c). For b = 6.0 
and 7.2 nm, spectra for a 1 meV mismatch are also given. 
One-to-one comparison between theory and experiment 
of absolute values is not helpful, because the QM's (QD's) 
actually behave in a very complex way [ p^ . In particular, 
25 can vary from device-to-device, and probably it de- 
creases with A^ [|l^. Nonetheless, the overall agreement 
between theory and experiment of the general spectral 
shape is quite good, indicating the crucial role played by 
mismatch. In particular, the appearance of the spectra 
in the weak coupling limit for small A^ values is now cor- 
rectly given, as well as the evolution with b of the peak 
appearing at A^ = 8 for intermediate coupling. A com- 
parison between panels (a) and (c) of Fig. ^ reveals that 
for smaller values of b (< 4.8 nm), for a reasonable choice 
of parameters {to, 5), mismatch does not produce sizeable 
effects. The reason is that the electrons are still rather 
delocalized, and distributed fairly evenly between the two 
dots. Exceptions to this substantial delocalization may 
arise only when both the constituent single QD states are 
magic, as discussed below, at intermediate coupling. For 
larger interdot distances, mismatch induces electron lo- 
calization. The manner in which it happens is determined 
by the balance between interdot and intradot Coulomb 
repulsion, and by the degree of mismatch between the sp 
energy levels, and so is difficult to predict except in some 
trivial cases for certain model parameters {uj,5). For ex- 
ample, a large mismatch compared to huj will cause the 
QD of depth Vq — 5 to eventually 'go away empty'. 

Finally, still assuming perfect coherency, a deeper the- 
oretical understanding of heteronuclear QM dissociation 
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can now be gained from analysis of the evolution with b of 
the sp molecular wavefunctions. Thus, for each sp wave- 
function (pniair, z,d) = Uniair, z)e^'^^^x<j ^6 introducc a 
z-probability distribution function defined as 

V{z) = 2tt drr[u{r,z)]^ . (1) 
Jo 

Figure I shows V{z) for (a) iV = 6, (b) = 8, and (c) 
= 12 (deeper well always in the z > region), each 
for several values of b. States are labeled as a, ±7r, . . . 
depending on the / = 0, ±1, ±2, . . . sp angular momen- 
tum, and indicate the spins. In each sub-panel, the 
probability functions are plotted, ordered from bottom 
to top, according to the increasing energies of the or- 
bitals. For each b, the third component of the total spin 
and total orbital angular momentum of the ground state 
are also indicated by the standard spectroscopic notation 
2^-+i|L^| with E,n, A,.. . denoting \L^\ = 0,1,2,... We 
conclude that: i) QM's dissociate more easily for smaller 
values of b, if they yield magic number QD's, as is the case 
for iV = 12 ^ 6 6 for 6 = 4.8 nm (c) or TV = 4 ^ 2 -f 2 
(not shown), for example, ii) Particularly for intermedi- 
ate values of b, not all orbitals contribute equally to the 
QM bonding, i.e., the degree of hybridization is not the 
same for all QD sp orbitals. See for example the tt and 
a states in the b = 4.8 nm panel of (a), iii) At larger 
b, dissociation can lead to Hund's first rule like filling in 
one of the QD's and full shell filling in the other dot. See 
for example the b = 7.2 nm panel in (a) for N = 6, which 
dissociates into 2 -|- 4. The same happens for the A^ = 10 
QM, which dissociates into 4 -|- 6 (not shown). In other 
cases, dissociation leads to Hund's first rule like filling in 
each of the QD's, as shown in the b = 12 nm panel of 
(b) for A^ = 8, which breaks into 4 + 4. In close anal- 
ogy with natural molecules, atomic nuclei, or multiply 
charged simple metal clusters homo and heteronu- 
clear QM's choose preferred dissociation channels yield- 
ing the most stable QD configurations, iv) Some config- 
urations are extremely difficult to disentangle: even at 
very large b, there can still be orbitals contributing to 
the QM bonding. A good example of this is the N — 8 
QM for 6 = 12 nm (b). 
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FIG. 1. Schematic diagrams of (a) mesa containing two 

vertically coupled quantum dots and (c) double quantum well 
structure, and (b) scanning electron micrograph of a typical 
circular mesa. 
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FIG. 2. (a) Calculated A2(Ai')/A2(2) for homonuclear QM's with different interdot distances, b. Also shown is the calculated 
reference spectrum for a single QD. (b) Experimental QM addition energy spectra, A2(A'^)/A2_(2), for several interdot distances 
between 2.5 and 7.5 nm. Also shown is an experimental reference spectrum for a single QD ||l8|. (c) Same as panel (a) but for 
heteronuclear QM's obtained using a 25 = 2 meV mismatch (dotted lines for b = 6.0 and 7.2 nm are for 25 = 1 meV). In each 
panel the curves have been vertically offset so that at A'^ = 2 they are equally separated by 0.5 units for clarity. All traces in 
panels (a) and (c) except 3.6 and 6.0 nm: HiK) marks cases where we could clearly identify Hund's first rule like filling within 
single dot, or bonding or anti-bonding states (constituent dot states). 
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FIG. 3. Calculated probability distribution functions V{z) (arbitrary units) as a function of z for the heteronuclear A = 6, 
A = 8, and A = 12 QM's (a), (b), and (c) respectively, using a 25 = 2 meV mismatch. 
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